Antiferromagnetic Order of Repulsively Interacting Fermions on Optical lattices 
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The Neel state in fermionic mixtures of two pseudospin speeies in an optical lattice is analyzed at 
low temperatures. Experimentally it remains a challenge to demonstrate antiferromagnetic correla- 
tions in ultracold fermionic quantum gases. We find that, while in balanced systems the Neel order 
parameter can point in any spatial direction, in imbalanced mixtures antiferromagnetism is strictly 
perpendicular to the quantization axis (i.e., the z-axis). Since, experimentally, one always has to 
assume some minimal imbalance this should have important consequences for ongoing experiments. 

PACS numbers: 03.75.Hh, 03.75.Kk, 71.10.Fd 



INTRODUCTION 



Ultracold quantum gases on optical lattices provide 
an interesting experimental environment for testing el- 
ementary quantum many-body models like the Bose- or 
the Fermi-Hubbard model [3, 0). In fact, microscopic 
quantum phenomena, predicted by theory, like correlated 
particle tunneling Q and superexchange [H, as well as 
macroscopic quantum phenomena like the Mott metal- 
insulator transition Q have by now been observed ex- 
perimentally. 

One experimental challenge for the near future is the 
observation of the Neel state in an ultracold fermionic 
mixture 0, [H, Q . Fermi mixtures are realized with neu- 
tral fermionic atoms in different hypcrfine states, for ex- 
ample 40 K with total angular momentum F = 9/2 in 
the hyperfme states F m — —9/2 and F m — —7/2, where 
the quantization axis is parallel to the magnetic field su- 
perimposed on the experimental environment [Io| . It is 
by now experimentally feasible to detect the two differ- 
ent hyperfme states seperately with in-situ imaging tech- 
niqu es [9] , [ill, 
ods lid, LL 




X 31 ] as well as with timc-of-flight meth- 



so that it seems a priori possible that the 
detection of antiferromagnetic correlations is observable 
with one of those pscudospin-selective detection meth- 
ods [15[ . Clearly a thorough understanding of theoretical 
concepts and predictions will be of great help in planning 
any experimental investigation of the antiferromagnetic 
state. 

In this paper we focus on the spatial structure of the 
Neel state in trapped ultracold Fermi mixtures. We show 
that antiferromagnetic correlations of a two spin-species 
Fermi mixture in the Neel state are expected to be spa- 
tially perpendicular to the z-axis (of the underlying Hub- 
bard model) in general, so that pscudospin-selective mea- 
surement methods may fail to detect such transverse an- 
tiferromagnetism. In this paper we first motivate our 
methods, then present our results and close with a pro- 
posal on how the transverse Neel state may be detected 
experimentally. 



MODEL AND METHOD 

A two pseudospin-species Fermi mixture on an optical 
lattice may be described by an inhomogeneous Hubbard 
model 

0,1 

TtI [l8j]. Since, for the purposes of this paper, 
we are interested in the phase diagram of the Hubbard 
model at low temperatures and not-too-strong repulsive 
interaction (0 < U < bandwidth rj At), it suffices to 
study the Hubbard Hamiltonian in the saddle point ap- 
proximation: 



(ij),<r io ~ 
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(1) 



where c icr is the fermionic annihilation operator for a 
fermion at site i with pseudospin a (e.g., f = F m = —9/2 
and | = F m = -7/2), n ia = c\ a c ii7 , n, = \{h- A + 

and S; is the local spin operator. Furthermore, t is 
the nearest-neighbor hopping amplitude, V is the con- 
fining potential strength and U describes the on-site re- 
pulsion. The spin-dependent chemical potential fi a con- 
trols the particle numbers in the grand-canonical formal- 
ism. The saddle point approximation ([1]) of the Hub- 
bard Hamiltonian is able to describe antiferromagnetic 
correlations both in z-direction and in the xy-plane. It 
is expected to yield qualitatively correct results at not- 
too-large coupling (U < bandwidth) but is known to 
quantitatively overestimate energy scales (gaps, critical 
temperatures) [l9j]. However, it is difficult to determine 
the accuracy of the saddle-point approximation exactly 
for a trapped inhomogeneous systems. Since the filling 
is space-dependent and the system has a finite size one 
would not expect the logarithmic van Hove singularity 
in the non-interacting density of states in two dimen- 
sions to cause a divergent contribution in the second or- 
der perturbation expansion [ljjj]. Hence, the results for 
two dimensions are expected to be similar to the ones 
in three dimensions and therefore it is our best guess is 
that the order parameter in the saddle-point approxima- 
tion is suppressed by a factor of two to four by quantum 
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FIG. 1: Particle density n (blue/dark line) and magnetizations 5*3 (orange/grey line) and Si (green/light grey line) on a 42 x 42 
square lattice depending on ^-position. Magnetic ordering perpendicular to the z-axis is suppressed in the upper panels while 
it is allowed in the lower panels. The parameters are chosen to be V = 0.01, U = 2.4, /3 = 10 and = 1.5 in units of t. 

The unbalance parameter A/x = /x-f — /xj. is chosen to be (a) A/x = , (b) A/x = 0.2 and (c) A/t = 0.6. 



fluctuations. 

For balanced systems, A/x = /Xj — /t j = 0, the Hamil- 
tonian (JTJ) is fully invariant under spin-rotations of the 
form 

CiT ) -> U f CiT ) , hi G 517(2) , (2) 

while the symmetry group reduces to U(l) for any A/x ^ 
0. Since one expects the grand potential of any self- 
consistent antifcrromagnctic solution to have the same 
global symmetry as the Hamiltonian (Q]), the Neel state 
should break a SU (2)-symmctry for A/x = and an U(l)- 
symmetry for A/x ^ 0. 

The Hamiltonian ((T|) can conveniently be written in 
the following matrix form: 

(3) 

where we defined S± = Si ± iS 2 , and TCoa contains the 
spin-diagonal terms of the Hamiltonian |T]). The self- 
consistency equations take the form: 

(K) = El«u,«.| 2 //9(^) (4) 

s=l 

(S i+ ) = ^2ut t1 u is J p (E s ) , (5) 

s=l 



where Ui SlCT are the components of an unitary transfor- 
mation diagonalizing ([3]), _E S are the eigenvalues of the 
saddle point Hamiltonian JT]), M is the total number of 
lattice sites and fp(x) = [1 + cxp(/3x)] _1 is the Fermi 
function at inverse temperature (5 = 1/fcsT. 

We determine antiferromagnetic solutions self-consis- 
tently by iterating the self-consistency equations ([3]) and 
(|5|) in combination with the Hamiltonian j3]). In order 
to find solutions with restricted antiferromagnetic order 
(purely in z-dircction) wc start the iteration cycle with 
(1/2 i) = ^ ^ or au ia ttice sites i in the initial state. This 
automatically brings the Hamiltonian ([3]) into a block- 
diagonal structure and, therefore, no antifcrromagnetism 
in xy-direction occurs at any iteration step, since in this 
case at least one of the coefficients u- ls ^ or Ui s ,x van- 
ishes for all indices {i, s}. We also solve the more general 
problem with possibly also (S x , 2 .) ^ and compare the 
grand potentials of both solutions at the same values of 
the parameters V, t, U, /Xj, /Xj and (3. 

RESULTS 

We now present a selection of numerical results for 
different spin-averaged chemical potentials fi = gC^T + 
Hi ) and different imbalance strengths on a 42 x 42 square 
lattice at fixed parameters V = 0.01, U = 2.4 and (3 = 10 
in units of the nearest-neighbor hopping amplitude t. In 
Figure [T] we choose /x = 1.5 so that an approximately 
half-filled region (i.e., with one particle per site) lies in 



3 






(a) 



(b) 



(c) 



FIG. 2: Particle density n (blue/dark line) and magnetizations 5*3 (orange/grey line) and Si (green/light grey line) on a 42 x 42 
square lattice depending on ^-position. Magnetic ordering perpendicular to the z-axis is suppressed in the upper panels while 
it is allowed in the lower panels. The parameters are chosen to be V = 0.01, U = 2.4, /3 = 10 and = 2.0 in units of t. 

The unbalance parameter A/i = /ij — /i; is chosen to be (a) A[i = , (b) A/z = 0.2 and (c) Afj, = 0.6. 



the center of the trap. In Figure [5] we increase the filling 
and choose [i = 2.0, so that the center of the trap is more 
than half-filled. Our results also confirm the suggestion of 
17 , HI , that antiferromagnetic correlations are enhanced 



in regions which are nearly half-filled. 

In general, we may obtain solutions with antiferromag- 
netic correlations in z-direction as well as in the xy-planc. 
Results for a system with a balanced number of particles 
are shown in Figures QJa) and^a). The spin-rotational 
symmetry of the Hamiltonian ([I]) at = fii is not bro- 
ken, and, accordingly, the amplitudes of the magnetiza- 
tion components (S 3 ) (upper panel, assuming (S^) = 0) 
and (Sj) (lower panel, choosing {S 3 ) = 0) are identical 
in Figures [Ha) and[DJa). The grand potentials of both 
solutions are, of course, degenerate. 

In unbalanced systems (/i| ^ /ij) the spin-rotational 
symmetry is reduced from SU(2) to U(l) and the Hamil- 
tonian |1} is invariant only under rotations around the 
z-axis. Results with the magnetization restricted to the 
z-direction are shown in the upper panels of Figures [1] and 
[21 and results allowing also for (S 1 ) ^ are shown in the 
lower panels. As one can see from Figures [T] and [51 anti- 
ferromagnetic order decreases as the imbalance (i.e., A/i) 
is increased. In Figures QJb) and [2f b) numerically stable 
solutions are presented with antiferromagnetic order in 
the z- and the x-direction, respectively. Note that the 
solution with its order parameter perpendicular to the 
z-axis has a lower grand potential and is therefore the 
thermodynamically stable one. In Figures [ljc) and[^c) 
results with a still stronger imbalance arc shown. While 



antifcrromagnctism purely in z-dircction is not numeri- 
cally stable in this case, staggered order still persists in 
the icy-plane and lowers the grand potential. 

We exclusively found solutions with S+j exp (—i<t>) £ R 
and a site- independent phase <p, so that the antiferromag- 
netic order parameter in the xy-plane can by a canoni- 
cal transformation always be chosen to be perpendicular 
to the y-axis. This result is a priori not at all obvious, 
but none of our numerical investigations yielded a self- 
consistent solution with a more complicated structure in 
the xy-plane. In the following we assume that the spins 
have been globally rotated such that 6*24 = for all i. 

Incommensurate states in the translationally invariant 
Hubbard model have been postulated in the literature 
already long ago [2(j ■ Here we find a somewhat different 
result, with a checkerboard structure for the and 
componcnts of the antiferromagnetic order parameter in 
several spatially distinct regions, which are separated by 
domain walls. These domain walls occur in regions where 
the total filling changes from ??i = (hi) « 0.5 to smaller or 
bigger values. We illustrate this domain wall formation in 
FigurcEl where we plot S14 cxp(zQ -i), with Q = (vr, 7r) T , 
as a function of the spatial coordinates x and y. A per- 
fectly commensurate Neel state in the balanced homoge- 
neous Hubbard model would not have any sign change. 
In contrast, we find sign changes in every thermodynam- 
ically stable self-consistent antiferromagnetic pattern in 
either balanced (A/i = 0) or unbalanced (A// 7^ 0) sys- 
tems. Furthermore, we find that any solution minimizing 
the grand potential systematically has a antisymmctrical 
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spatial structure in quantities describing antifcrromag- 
netic correlations (e.g., Sij at any A/i- value or S^i at 
A/z = 0), while the other quantities (e.g., n\ and 5 3j j at 
Afi 7^ 0) have a spatially symmetrical structure. 



EXPERIMENTAL CONSEQUENCES 

We expect spin-selective measurement methods like in- 
situ imaging or time-of-flight (TOF) techniques to be in- 
adequate [15j for the successful detection of the Neel state 
in ultracold Fermi mixtures on two-dimensional optical 
lattices, since the antiferromagnetic order parameter is 
predicted to be perpendicular to the z-axis at any fi- 
nite imbalance. While the experimental techniques men- 
tioned above are able to distinguish between two hyper- 
fine states (e.g., F m = —9/2 and F m = —7/2), which 
corresponds to the measurement of S3 = \{n^ — njj, 
they are not able to capture coherent linear combina- 
tions represented by anomalous expectation values like 
Si = (S+ + S_)/2. 

Recently a method for reconstructing the non-diagonal 
elements of a real-space Green function through a combi- 
nation of a double pscudospin-mixingRaman pulse and 
a TOF measurement was proposed 21|. However, this 
proposal for detecting transverse components of the order 
parameter, too, seems problematic due to the remaining 
[/(l)-symmctry of the antiferromagnetic state. As a con- 
sequence of this symmetry, the direction of the transverse 
AFM order parameter (and, hence, the non-diagonal el- 
ements of the Green function) will be different in each 
TOF image, effectively preventing or at least seriously 
complicating the reconstruction of the Green function. 

We, therefore, propose to follow a different strategy: 
An observable that is readily accessible using in-situ 
imaging techniques @,[ll, 12, 13[ is the spatial distribu- 
tion of particle numbers n\. Our calculations show that 
antifcrromagnctism leads to a broadening of the approx- 
imately half-filled region (n\ ~ 0.5), even if the interac- 
tion is weak. This signature can be seen very clearly in 
Figures [T] and [2] which show that the antiferromagnetic 
states have a plateau in the half-filled region ( "wedding- 
cake structure"), that is characteristically absent in the 
paramagnetic state [upper plot of the (c)-series]. Our 
proposal is to use this plateau in the density profile as a 
smoking gun for antifcrromagnctism. In doing the experi- 
ment, it will of course be crucial to choose the parameters 
of the experimental setup such that the plateau becomes 
visible within the resolution of in-situ imaging. Since re- 
gions of sizes down to approximatively 5x5 lattice sites 
are resolvable, parameters similar to the ones chosen in 
Figure [T] seem to be adequate for this purpose. Obvi- 
ously, a plateau at half-filling might also arise at strong 
interaction (for U > U c > bandwidth) due to a Mott 
metal-insulator transition within a paramagnetic state 
H^. Hence, it would not be possible to distinguish 




FIG. 3: Staggered antiferromagnetic order parameter m = 
Si,i exp(iQ • i), with Q = (tv, 7r) T , on a 42 x 42 square lattice. 
The sign of m is positive (negative) in the dark (light) re- 
gions, indicating the formation of different antiferromagnetic 
domains, separated by walls of parallelly aligned spins. In 
particular, the sign of m is positive in the region of approx- 
imate half filling (ni ~ 0.5) and changes where the gradient 
of the filling factor is maximal. The parameters are chosen 
to be V = 0.01, U = 2.4, (3 — 10, A/j, = and = 2.0 

in units of t. The sign changes are also characteristical for 
unbalanced systems (Afj, 7^ 0). 



clearly between paramagnetic Mott-insulating states and 
antiferromagnetic ones in the strong-interaction regime. 
Our proposal is, therefore, to try to detect the broaden- 
ing of the approximately half-filled region as a signature 
for antiferromagnetic correlations in the weak- coupling 
regime (U < U c ), where the plateau cannot be caused by 
a Mott transition. This strategy should be practicable, 
since the interaction strength U in ultracold quantum 
gases is straightforwardly tunable experimentally Q ■ 



SUMMARY AND CONCLUSIONS 

We investigated a two pseudospin-species Hubbard- 
modcl well suitable to describe spatially inhomogeneous 
Neel states of Fermi mixtures on a two-dimensional op- 
tical lattice in a trap. In our saddle point approximation 
we found the antiferromagnetic order parameter to have 
an arbitrary direction in a balanced mixture [in accor- 
dance with the full SU (2)-symmetry of the Hamiltonian] 
and to be strictly perpendicular to the quantization axis 
at any imbalance strength [in accordance with the U(l)- 
symmetry for Afi 7^ 0]. Experimentally, one always has 
to assume some minimal imbalance. In order to exper- 
imentally observe these staggered magnetic structures, 
we propose to use in-situ imaging of the spatial distri- 
bution of particle numbers and to use "wedding-cake" 
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structures as a possible signature for antifcrromagnetism 
in the parameter regime corresponding to U < U c , where 
U c is the critical interaction for Mott-insulator transi- 
tions. Furthermore, we found that the antiferromagnetic 
order parameter may always be chosen to be parallel to 
a globally defined direction and that there is a forma- 
tion of domain-walls between checkerboard-type regions 
of the magnetization. The spins are aligned parallelly 
along these domain walls, which occur in regions where 
the filling factor changes from half-filling to higher or 
lower values. 

We thank Sebastian Will (Univ. Mainz) for interesting 
discussions on the possible experimental observation of 
antiferromagnetic order as proposed in this paper. 



[3] 



* Electronic address: tobias.gottwald@uni-mainz.de 
[1] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner and P. 

Zoller, Phys. Rev. Lett. 81, 3108 (1998). 
[2] I. Bloch, J. Dalibard and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

S. Foiling, S. Trotzky, P. Cheinet, M. Feld, R. Saers, 
A. Widera, T. Miiller, and I. Bloch, Nature 448, 1029 
(2007). 

[4] S. Trotzky, P. Cheinet, S. Foiling, M. Feld, U. Schnor- 

rberger, A. M. Rey, A. Polkovnikov, E. A. Demler, M. D. 

Lukin and I. Bloch, Science 319, 295 (2008). 
[5] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansen and 

I. Bloch, Nature 415, 39 (2002). 
[6] R. Jordens, N. Strohmaier, K. Giinter, H. Moritz and T. 

Esslinger, Nature 455, 204 (2008). 



[7] F. Werner, O. Parcollet, A. Georges and S. R. Hassan, 

Phys. Rev. Lett. 95, 056401 (2005). 
[8] A. Koetsier, R. A. Duine, I. Bloch and H. T. C. Stoof, 

Phys. Rev. A 77, 023623 (2008). 
[9] U. Schneider, L. Hackermuller, S. Will, T. Best, I. Bloch, 
T. A. Costi, R. W. Helmes, D. Rasch and A. Rosch, Sci- 
ence 322, 1520 (2008). 
[10] M. W. Zwierlein, A. Schirotzek, C. H. Schunck and W. 

Ketterle, Science 311, 492 (2006). 
[11] G. B. Partridge, Wenhui Li, R. I. Kamar, Yean-an Liao 

and R. G. Hulet, Science 311, 503 (2006). 
[12] G. B. Partridge, Wenhui Li, Y. A. Liao, R. G. Hulet, M. 
Haque and H. T. C. Stoof, Phys. Rev. Lett. 97, 190407 
(2006). 

[13] K. B. Gubbels, M. W. J. Romans and H. T. C. Stoof, 

Phys. Rev. Lett. 97, 210402 (2006). 
[14] F. Gerbier, A. Widera, S. Foiling, O. Mandel, T. Gericke 

and I. Bloch, Phys. Rev. A 72, 053606 (2005). 
[15] G. M. Bruun, B. M. Andersen, E. Demler and A. S. 

S0rensen, Phys. Rev. Lett. 102, 030401 (2009). 
[16] M. Iskin and C. J. Williams, Phys. Rev. A 78, 011603(R) 

(2008). 

[17] B. M. Andersen and G. M. Bruun, Phys. Rev. A 76, 

041602(R) (2007). 
[18] M. Snoek, I. Titvinidze, C. Toke, K. Byczuk and W. 

Hofstetter, New J. Phys. 10, 093008 (2008). 
[19] T. Schauerte and P. G. J. van Dongen, Phys. Rev. B 65, 

081105(R) (2002). 
[20] H. J. Schulz, Phys. Rev. Lett. 64, 1445 (1990). 
[21] L.-M. Duan, Phys. Rev. Lett. 96, 103201 (2006). 
[22] S. Foiling, A. Widera, T. Miiller, F. Gerbier and I. Bloch, 

Phys. Rev. Lett. 97, 060403 (2006). 
[23] R. W. Helmes, T. A. Costi and A. Rosch, Phys. Rev. 

Lett. 100, 056403 (2008). 



